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Abstract 

The decay rates of metastable SQCD vacua in ISS-type models, both towards su- 
persymmetric vacua as well as towards other nonsupersymmetric configurations 
arising in theories with elementary spectators, are estimated numerically in the 
semiclassical approximation by computing the corresponding multifield bounce 
configurations. The scahng of the bounce action with respect to the most relevant 
dimensionless couplings and ratios of scales is analyzed. In the case of the decays 
towards the susy vacua generated by nonperturbative efi^ects, the results confirm 
previous analytical estimations of this scaling, obtained by assuming a triangular 
potential barrier. The decay rates towards susy vacua generated by R-symmetry 
breaking interactions turn out to be more than sufficiently suppressed for the 
phenomenologically relevant parameter range, and their behavior in this regime 
differs from analytic estimations valid for parametrically small scale ratios. It 
is also shown that in models with spectator fields, even though the decays to- 
wards vacua involving nonzero spectator VEVs don't have a strong parametric 
dependence on the scale ratios, the ISS vacuum can still be made long-lived in 
the presence of R-symmetry breaking interactions. 



1 Introduction 



Since Intriligator, Seiberg and Shih showed in their seminal paper [T] that long lived, 
supersymmetry-breaking metastable vacua arose in simple theories, such as SQCD in the 
free magnetic phase, many works have studied models of particle physics based upon ISS- 
type constructions. A minimum requirement for these models to be considered as plausible 
descriptions of the particle interactions in Nature is that the lifetime of the relevant vacuum 
should be greater than the estimated age of the Universe. 

Ref. argued that the lifetime of the ISS vacuum in the original model depended 
parametrically on the ratio between the dynamical scale A of the gauge interactions and 
the mass scale n of the low-energy magnetic theory, which is related to the masses of the 
electric quarks. Thus, long enough lifetimes are guaranteed for fi/A ^ 1, which can always 
be achieved by pushing A to higher energies. A similar situation holds for models in which R- 
symmetry violating interactions involving a new scale /x,^ are introduced, as in ref. in this 
case the decay towards supersymmetric vacua is parametrically suppressed for /i^Z/x -C 1. 
However, since //^ determines the gaugino masses while fi ends up controlling the squark 
soft masses, the ratio /i,/,//^ cannot be made too small without creating a very fine-tuned 
spectrum. In principle it is not completely clear whether this will be in conflict with the 
vacuum being sufficiently stable, but the results of ref. showed that this is not necessarily 
the case. Here a more a more detailed numerical analysis will be presented. 

Another threat to the stability of the vacuum appears when the ISS model is coupled 
to elementary spectators. In models based upon ISS constructions and attempting to link 
supersymmetry breaking with flavor physics, such as those of refs. O IH [5l [6], some of the 
generations of the Standard Model are embedded into the meson fields of the low energy 
theory; generically, this cannot be done without having light mesons in undesired representa- 
tions of the gauge group, which have to be lifted by coupling them to elementary spectators. 
These new couplings give rise to new metastable vacua involving nonzero expectation values 
for the spectator fields, which have lower energy than the ISS vacuum]^ In contrast 
with the decays towards supersymmetric minima, which exhibit the parametric suppres- 
sions alluded to before, it is not clear whether the decay rates towards the new vacua will 
have a strong dependence on the ratios of scales; at least for weak R-symmetry breaking 
interactions, the rates are expected to be mainly determined by a dimensionless coupling 
which cannot be made parametrically small in phenomenologically acceptable models. It 
was argued in refs. |5l [6] that by making the electric quarks appropriately non-degenerate, 
the energy of the new vacua can be lifted above that of the ISS vacuum, which would pre- 
vent the decay. However, in models attempting to explain gaugino masses this solution is 
unsatisfactory, since the consideration of nondegenerate masses for the electric quarks leads 
to a suppression of the gaugino masses arising from R-symmetry breaking interactions [6]. 

^These vacua were initially found by D. Green, A. Katz and Z. Komargodski. Spectator fields may also 
be added to avoid tachyonic squarks in dynamical realizations of minimal gauge mediation, as in ref. [7], 
though in this case they do not give rise to lower energy vacua. Other models use spectator fields to build 
dynamical models of gaugino mediation [8|. 
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Thus, in models with degenerate quark masses and spectator fields, with no extra mass 
scales involved in the interactions of the latter -otherwise it is straightforward to lift the 
vacua involving spectator VEVs above the ISS vacuum- a more detailed estimation of the 
decay rates is needed in order to discern whether the ISS-type vacua are long-lived enough. 

This paper presents new calculations of the decay rates of ISS vacua in models with 
degenerate electric quarks, both towards supersymmetric vacua and towards vacua involving 
expectation values for spectator fields, in theories with and without R-symmetry breaking 
interactions. These calculations were performed by solving numerically for the semiclassical 
bounce configurations that dominate the path integral representing the decay amplitude [U] . 
The computation is not straightforward due to the fact that the bounce configurations involve 
several fields; it was tackled following the approach of ref. [10], with up to five fields. Also, 
since some of the vacua are only stabilized after taking into account quantum corrections, 
these have to be taken into account when solving for the bounce configurations. Partial 
results were already quoted in ref. [H] ; here a more complete presentation is given as well as 
results involving less approximations. 

The results confirm the qualitative dependence of the bounce action on dimensionless 
couplings and ratios of scales that can be extracted by identifying the small parameters that 
dominate the relevant contributions to the potential. Thus, the decays towards susy vacua 
generated by nonperturbative interactions are indeed parametrically suppressed for growing 
values of the ratio ^. Similarly, the decays towards susy vacua generated by R-symmetry 
breaking interactions are suppressed for diminishing Finally, the decays towards vacua 
involving nonzero VEVs of spectator fields are controlled essentially by the dimensionless 
parameter h present in the ISS model superpotential, although for fixed h they are enhanced 
by increasing This last property could not be captured by the simplified calculations of 
ref. in], which ignored R-symmetry breaking interactions, and used only a rough estimation 
of the gauge fields' one-loop contributions to the scalar potential. Analytic estimations of 
the scaling of the bounce action with respect to the above dimensionless parameters, that 
were previously obtained in the literature by assuming a triangular barrier and restricting 
to one field direction, are confirmed for the decays towards the nonperturbative susy vacua. 
In the case of susy vacua with fifp 0, the analytical scaling obtained in ref. [2j is valid 
for yU^Z/i <^ 1, which is not the regime studied here -numerical fiuctuations make it hard 
to explore. The scaling obtained here differs from the triangular barrier result, but the 
latter could still be valid for fi^/ fi <C 1. Regarding the resulting lifetimes, all decays to susy 
vacua, including the case with R-symmetry breaking interactions, turn out to be more than 
sufficiently suppressed for typical, phenomenologically acceptable values of the parameters, so 
that demanding long-lived vacua does not impose phenomenologically relevant constraints. 
In the case of models with spectator fields, centering on the case in which the magnetic 
group is trivial and for diagonal spectator fields (i.e. coupling to the diagonal of the meson 
field matrix), sufficiently long lifetimes can be achieved for small, nonzero fifp (perfectly 
safe regarding the decays towards susy vacua) and order one values of h; for the choice of 



parameters of [3.5, the bounce action for sufficiently high values of ^ > 4 ■ 10 ^ scales 
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roughly (up to errors around 10%) as 



clSS^spect -J P0 /-I -I X 

Regarding off-diagonal spectators (coupling to off-diagonal mesons), they are expected to 
generate milder instabilities than the diagonal ones: the new vacua are either further from 
the origin of field space or involve more fields getting a VEV, which implies an increase 
in the number of tunneling directions. Hence, if a model is safe with respect to decays 
towards vacua associated with diagonal spectator fields, it is expected that it will also be 
safe concerning decays to vacua generated by off-diagonal spectators. Thus only the first 
type decays were analyzed in detail. 

Finally, it should be commented that there exist other potential sources of instabilities 
in ISS-type models, which do not involve tunneling processes but rather growing classical 
solutions, such as Q-balls. Ref. [TT] argued that by making the dimensionless coupling h 
of the ISS model greater than the gauge coupling g of the gauged global symmetry, h > g, 
then Q-balls are unstable towards decay to massive vector bosons and they do not threat 
the stability of the theory. This requirement seems in tension with the fact that the decay 



rate towards spectator vacua grows for increasing values of h -see eq. 1.1 The above bound 
applies when the full global symmetry is gauged, so that the Goldstone modes of the global 
symmetry are eaten away and do not support Q-ball solutions. The situation worsens in 
models in which only a subgroup is gauged, since then the Goldstone modes become pseudo- 
Goldstone modes with masses of the order g/4:7c, which can support Q-balls and cannot decay 
to other particles. However, there is a possible way out consisting in gauging the global U(l) 
symmetry that guarantees the existence of Q-balls. This gauging seems to make the Q-balls 
unstable, at least for large values of the gauge coupling [121 113]. 

The organization of the rest of the paper is as follows. Section |2] reviews the ISS model, 
including the ISS vacuum, the supersymmetric vacua and the vacua arising when spectator 



fields are added. Section |3] deals with the numerical calculations. ^3J^ introduces the numer- 



ical method, while ^3.2 introduces a dimensionless parameterization of the bounce action, 
which allows to extract qualitative conclusions about the dependence of the latter with re- 
spect to the couplings and scale ratios of the theory. The results for decays towards susy 
vacua with and without R-symmetry breaking and for vacua involving spectator fields are 



given in §^3.4, 3.3 and 3.5, respectively. An appendix has been included, providing formulae 



for the tree-level mass matrices and the Coleman- Weinberg potential. 



2 Vacuum configurations in ISS-type models 

This paper focuses on weakly coupled ISS-type supersymmetric models involving an SU{Nc) 
gauge group, an Nf x Nj matrix $ of singlet meson fields, and fields q and q -"magnetic 
quarks"- in the fundamental and antifundamental representations of the gauge group. The 



3 



superpotential for the basic model is P 

W^ag = -Mr(/i2$) + /itr(g$g), (2.1) 

where the suffix "mag" is chosen to emphasize the fact that these models arise as low energy 
descriptions of SU (N) SQCD with Nj massive flavors in the free magnetic phase, Nc < 
Nf < |iVc, with Nc = Nf — Nc [H]. /i is a matrix of mass parameters, which can be related 
to the mass matrix of the electric quarks in the dual theory as 

—hfi'^ ~ Am, 

A being the dynamical scale of the electric SU (N) gauge group. For the rest of the paper, 
the focus will be on degenerate electric quarks, yielding fi = /iIat^. 

Apart from the weakly coupled superpotential Wmag, one has to take into account non- 
perturbative contributions generated by the gauge dynamics, 

Wnp = Ncih^^'AT'-''' det (2.2) 
where A^ is the dynamical scale in the magnetic theory, related to A as [T] 

The theory has a global SU {Nf) x SU {Nf) x U {l)v symmetry under which the fields trans- 
form as 





SU{Nf) 


SU{Nf) 


U{l)v 
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0. 



In typical models, as in applications to gauge mediation or models with composite gener- 
ations, a diagonal subgroup SU{Nk)D, < Nf is gauged, where SU{Nk) contains the 
Standard Model group. Similarly, models usually involve the gauging of one or more U{1) 
factors in order to lift massless Goldstone modes around the soon to be reviewed ISS vacuum 

i- 

Furthermore, the possibility of coupling the meson field to a matrix S of elementary spec- 
tator fields will be considered, as in the models of refs. [31 IH El E], without the introduction 
of further mass scales. This entails a superpotential contribution of the form 

Inspect = AATr$S. (2.3) 

To avoid spoiling the ISS vacuum, in which $ acquires nonzero F-terms along the diagonal, 
S is taken as a traceless matrix. If q and q are chosen to be in a reducible representation 
of SU{Nk)D, S can have components transforming in both vector representations of this 
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group -along diagonal blocks- and chiral ones -along off-diagonal blocks. Finally, in models 
of supersymmetry breaking based upon ISS vacua and aiming to explain gaugino masses, it 
is necessary to introduce an explicit source of violation of R-symmetry. We will follow the 
approach of ref. ^ and add the following contribution to the superpotential: 

WR = ^h^fi^Tr^\ (2.4) 

Throughout the rest of the paper it will assumed that the Kahler potential for the fields 
is canonical; this is justified since the models considered arise as weakly coupled infrared 
descriptions of the strongly coupled electric theories. 



2.1 Vacua in the R-symmetric case 

This section provides a brief review of the vacuum configurations of models without R- 
symmetry breaking, with superpotential of the form 



W = W^^g + Wnp + Inspect, 



(2.5) 



where the different terms are given in eqs. (2.1), (2.2) and (2.3). 



2.1.1 Susy vacua 



The superpotential of eq. (2.5) admits supersymmetric configurations, thanks to the non- 
perturbative contributions. These are jl] 



3Nc 



V 

A* \ Nf-Nc 



.A, 

with the rest of the fields having zero VEVs. 



2.1.2 ISS vacuum 



Ignoring the nonperturbative contributions Wnp to the superpotential, which can be done 
when the fields take values near the origin, the ISS model breaks supersymmetry due to the 
rank condition: 

- FI = -hfxH + hqq^ y^O, (2.6) 

because the two terms in the sum do not have the same rank; at most one cancel rank(gg) = 
Nc < Nf F-terms. Thus, supersymmetry is broken. Decomposing the meson and magnetic 
quarks in the following way. 
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q = I '^NcxNc j q = ( ^^cxffc 

yP{Nf-ffc)xffcJ ' \P{Nf-Nc)xNc 

then the resulting ISS vacuum is characterized by nonzero F-terms Fx = — /i/x^, with 

{XX)=P%^, (2.7) 
and the rest of the fields having zero VEVs. The tree-level vacuum energy is 

To lowest order, X is a pseudomodulus, but it is stabilized at the origin by one-loop quantum 
corrections It should be emphasized that all field directions, including the Goldstone 



modes associated with the global symmetries broken by the configuration of eq. (2.7), are 
stabilized once one gauges the global U(l) symmetry of the model [2]. Thus, the decays to 
other vacua can only proceed through tunneling (as long as the Q-balls of ref. [11] are made 
unstable, as commented in the introduction). 



2.1.3 Vacua involving spectator VEVs 

As was motivated before, the spectator fields can be embedded into a traceless matrix. They 



are used to lift some mesons included in the field X of eq. (2.1.2), which are light in the 



ISS vacuum; thus, 5* only couples to X. Taking into account the spectator couplings, and 



ignoring again nonperturbative contributions to the superpotential, instead of eq. (2.6) one 
has 

- fI = -hfx%^ + hxx\ -Fj^ = -hixHN^ + hpp^ + XAS. (2.8) 

The rank condition applies as before to the fields X but not to p, p; one can use the 
spectators to cancel more F-terms than in the ISS vacuum. [5]. 

One can then consider VEVs for both off-diagonal as well as diagonal components of the 
spectator fields. In the case of diagonal VEVs, if S is an Nj x Nj traceless matrix coupling 
to an Nj X Nj block Xj of the meson field, one can use Nj — 1 diagonal elements to cancel as 
many F-terms of Xj, and then the uncancelled F-term can be made zero by giving a VEV 
to a pair of components of the p, p fields. This in turn requires to set one component of X 
to zero in order not to turn on some components of Fz, F^- The new vacuum is [5] 

ixx) = {pp)jxj = /^'diag(0, . . . , 0, j), (S) = ^diag(l, . . . , 1, -(j - 1)), (2.9) 

whose tree-level energy is 

Aspect = /i'(iVc-J + l)/i', (2.10) 

i.e., lower than the energy of the ISS vacuum. In this vacuum, Y as well as the components 
of X that do not couple to spectators are tree-level fiat directions; again, they are stabilized 
by quantum corrections. 
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Regarding vacua with VEVs for off-diagonal spectators, one may distinguish two cases. 
The first case is that in which off-diagonal spectator fields come in conjugate pairs. This 
allows again to cancel more F-terms than the Nc ones allowed by the rank condition for the 
ISS vacuum. Considering p, p a.s Nc x Nc matrices, one can give identical VEVs to multiple 
rows of p,p; this allows to cancel the diagonal F-terms in eq. (2.8) but generates symmetric 
off-diagonal contributions, which can be compensated in turn by giving VEVs to off-diagonal 
spectators. The F-term conditions for Fz, then force to set to zero one linear combination 
of the rows of the Nc x Nc matrix of VEVs of x and x- Hence, if one cancels j F-terms along 
X, one needs to give VEVs to j rows of p, p as well as ^^"^j^ pairs of off-diagonal spectators; 
the resulting vacuum energy is again that of eq. (2.10). 

Alternatively, one can have off-diagonal spectators that are not accompanied by their 
conjugates. Again, one can cancel the diagonal F-terms of X as in the previous paragraph 
by giving VEVs to the p, p fields, which generates the same pairs of off-diagonal contributions 
to Fx- For each of these pairs, one can cancel one of the F-terms by giving a VEV to an 
off-diagonal spectator field; since it is assumed that the spectators do not come in conjugate 
pairs, there will remain a nonzero off-diagonal contribution, — -^xij — h Pip] ^ 0, i ^ j. As 
shown in [6] , the resulting vacuum configurations have a runaway towards vacua with energy 
as in eq. (2.10). To see this, consider the cancellation of a pair of diagonal components of 
-FI, say -Fj^.. and -Pj^.., by choosing 

PipJ = P,p] = p''. (2.11) 

An off-diagonal spectator Sji can be given a VEV to cancel then {Fj^)ij = hpiPj is 

in principle nonzero but can be minimized by taking either pi or pj to zero, which by virtue 
of (2.11) implies that either pi or pj must approach infinity, which generates a runaway. 

At this point it should be clear that the tunneling from the ISS vacuum towards the 
vacua of this section should proceed preferably towards those involving spectators along the 
diagonal, either because they are closer to the ISS vacuum or because the decay involves a 
smaller number of tunneling directions. For this reason the numerical estimates in ^ center 
on the case of diagonal spectators. 



2.2 



Vacua in the presence of R-symmetry violating interactions 



This section is dedicated to the vacuum configurations in the presence of R-symmetry break- 



ing interactions coming from the superpotential contribution of eq. (2.4), so that the total 
superpotential is of the form 



np 



where the different terms are given in eqs. (2.1), (2.2), (2.3) and (2.4) 



In this case, the violation of R-symmetry gives rise to susy-vacua near the origin, as 
expected from the results of ref. [I5]. These vacua are given by 
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If the R-symmetry breaking deformation W^, of the superpotential is small, it is expected 



that there will be vacuum configurations deforming the ISS vacuum of §2.1. 2[ This turns 
out to be the case indeed; although at tree-level the pseudomodulus field X of the ISS model 
becomes unstable, one-loop quantum corrections stabilize it at a nonzero VEV, so that the 
new vacuum becomes [2j 



' c 



Analogously, the vacua involving spectator fields also survive the /i,^ deformation; in this 
case the pseudomodulus field Y gets a nonzero VEV due to quantum corrections, of the same 



order as X in eq. (2.12). More details about the Coleman- Weinberg potential responsible 



for these quantum corrections are given in 



3 Numerical estimations of decay rates 
3.1 Method 

In order to estimate the decay rates of the ISS vacuum towards the supersymmetric and 
nonsupersymmetric vacua describe above, the semiclassical approach of Coleman [H] will 
be followed. In this approach, the lifetime of a metastable vacuum is computed from a 
saddlepoint evaluation of the path integral for the transition amplitude of the false vacuum 
onto itself after infinite time. In a theory involving a set of fields (p = {(pi}, the decay rate 
per unite volume can be written as [9l [16] 

-1/2 

x{l + 0{h)), 5fe = 5(0)-5(0_). (3.1) 



r _s. Si det'[-D + V"(''' 

— = p h — 

V 4772^2 



det[-n + K"((/>_)] 

In the previous formula, 0_ denotes the metastable vacuum configuration, while det' indicates 
the determinant with translational zero modes omitted. represents the semiclassical bounce 
configuration, which is obtained as a solution of the Euclidean equations of motion with 
appropriate boundary conditions. With the Euclidean action written as 



^[0.] = j d'x^- J2i9,<l>^y + V^(0.), (3.2) 



then the bounce configuration satisfies the equations 

dV((f)) 

0(Pi T->00 

Assuming an 0(4) symmetry, writing r = (r, x), the equations above turn into 
c?20i 3d^i dV{(j)) , , 

-— + -— = ——, lim0 = 0_, —0 = 0. 3.3 
ar-^ r dr ochj r-s>oo dr 
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In the case of a single field, eq. (3.3) can be solved numerically by employing a shooting 
method. However, this is not suited for multifield configurations. The calculations pre- 
sented in this paper were done with the technique of ref. [10]. Essentially, this approach 



contemplates the equation (3.3) with a modified friction term 



<f(Pi a-ld(f)i dV{4>) 



+ 



dr 



(3.4) 



The equations for the frictionless case a = 1 are equivalent to the equations of motion of the 
following one-dimensional action 



S 



dr 



(3.5) 



However, the bounce solution is not strictly a minimum of the above functional, but rather 
corresponds to a saddlepoint |IB] -solutions staying closer to the true vacuum than the 
bounce have lower values of 5*. Nevertheless, by using a modified potential K, obtained 
from V by flattening out the potential energies below that of the false vacuum, it can be 
argued that the frictionless bounce configuration for a = 1 (see eq. (3.4)) is a minimum of 
the functional obtained by substituting V with in eq. (3.5), in the limit e — )■ [TU]. 
can be taken as 



r(0) - V(0_ 



(r(0) - v{<p.)y 



4 



+ e2 - 2e, 



+ 3e. 



where 0+ designates the lower lying vacuum. The reason for this is that, for e — > 0, the 
frictionless bounce has the same value of than field configurations joining both vacua 
and hence minimizing the functional; the bounce can be obtained from them by simply 
eliminating the points with V{(j)) < V{(f)^), whose contribution to the action vanishes for 
e ^ 0. 

Once the solution of eq. (3.4) for a = 1 has been obtained via a minimization procedure, 
the sought for solution corresponding to a = 4 can be obtained iteratively in incremental 



steps in a by solving a finite-difference version of eq. (3.3) with a linearized gradient of the 
potential. Discretizing the r variable as an n point lattice = kAr, fc = 0,...,n — 1, and 
denoting (pkifi) = 4>k^ linearized equation for a given value of a and a previous estimation 
(/) of the solution is 



Ar2 



+ 



a — 1 



'i + Az)Ar 2Ar 



d^yi^') + (0i - ~^i)didkV{^' 



which can be easily inverted. Ai is an offset parameter introduced to avoid pathological 
behavior near z = 0, which eventually has to be driven to zero. In the computations presented 
here, the iterations were performed with increasing number of steps in a until the successive 
values of the bounce action differed in less than a few percent. 
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In the previous formulae, V{(f)) is in principle meant to be the tree-level potential. How- 
ever, when the higher energy vacuum is only stabilized after taking into account quantum 
corrections -as happens in the /i,^ deformation of the ISS vacuum- there will be no bounce 



solution to eq. (3.3). An obvious way out is to substitute the tree-level potential by the 
perturbative quantum effective potential Kff(</') [IZj- This approach is correct as long as the 
effect of higher derivatives terms or field renormalizations can be ignored; to take these into 
account, and to avoid problems with imaginary values of the effective potential, a modified 
formalism is needed, such as the one described in ref . [18] . To leading order both approaches 
are expected to be equivalent; therefore in the computations presented here it was pro- 
ceeded by simply using the one-loop effective potential, which was obtained by evaluating 
the Coleman- Weinberg potential for arbitrary background values of the fields, as detailed in 

1 

The numerical method summarized above allows to compute the bounce action S{(f)) in 



eq. (3.1). The calculation of the prefactor is a much more daunting task; however, since 
the bounce action enters exponentially, its contribution is dominant. The prefactor can be 
estimated as an appropriate power of the most relevant scale in the potential; with scales 
of the order of 200 TeV -the typical mass scale for the fields p, Z in phenomenological 
models of direct mediation based on the ISS model, then a bounce action of the order 
Sh = S{(t>) — > 400 (in natural units) is enough to guarantee a lifetime of the order of 

the age of the Universe. 



3.2 Dimensionless parametrization of the bounce action 

Since the action of the bounce configuration is dimensionless in natural units, it depends 
only on dimensionless couplings or ratios of the scales in the potential. Thus, it is desirable 
to perform the computations using dimensionless variables. This not only might help avoid 
numerical fluctuations in the presence of very different mass scales, but also allows to reach 
qualitative conclusions about the dependence of the bounce action on the dimensionless 
quantities. 



The parameters in the superpotentials of eqs. (2.1), (2.2), (2.3) and (2.4) are the dimen- 



sionless couplings h, A and the scales /i,//^, A. From these we can define dimensionless scale 
ratios and fields, as well as a dimensionless radial variable and potential, as follows 



_ fi^ _ A 

/i /i 



= -0, f = ixr, 1/(0; h, A, Kx) = 



V{(j); /i,A,/X0,A). 



Then the "hatted" fields corresponding to the bounce configuration satisfy the following 
equation involving only dimensionless quantities, 



flf2 f. (If. 



lim ( 

r— >oo 



dr 



(0) = 0, 
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whereas the bounce action can be written as 



S = 2n' 



f df 



dr 



One can get some quahtative insight on the dependence of the bounce action on some of 
the dimensionless parameters when these get small [18]. Suppose that one can single out a 
coupling 5 such that V 
bounce action is 



\V is of order one. Rescaling the radial variable r = AjAf, the 



S 



27r2 



r'^dr 



dr 



+ V{<P) 



where the integrand is expected to be of order one, so that the bounce action will scale like 
1 

<5- 

In the case of decays from the ISS vacuum, generically h is assumed to be of order one; 

since it appears in the F-terms and the masses of the messenger fields, it cannot be made 

too small in phenomenologically acceptable models. The parameters which can be made 

small are k,^ and which are expected to determine the decay rates towards the susy 

vacua. This is so because both parameters appear in the superpotential contributions Wnp 

and Wr of eqs. (2.2) and (2.4), which are responsible for the existence of the susy vacua of 

1 



^2.1 



and ^2.2 



These superpotentials give rise to terms with positive powers of ^ and 
in the potential, and the discussion in the previous paragraph applies. Thus, the bounce 
action is expected to be parametrically large for small values of these ratios of scales, which 
guarantees a parametrically small decay rate. The same conclusions were achieved in refs. [1] 
and [2] using more detailed arguments based on one-dimensional tunneling with a triangular 
potential barrier, which yielded analytic estimates for the scaling. Ref. [1] argued that the 
bounce action associated with the decay towards the nonperturbative susy vacua should 
behave as 



A(Nf 



-3Nc) 

-Nc 



(3.6) 



while in ref. |2] the bounce action corresponding to the decay towards the susy vacua in the 
presence of R-symmetry breaking was estimated, in the limit <^ 1 as 



(3.7) 



The next section will present numerical results which confirm the qualitative behavior and 



the scaling of eq. (3.6). Regarding the decays to susy vacua in the presence of R-symmetry 
breaking, numerical fluctuations prevented to probe the regime <^ 1 for which eq. (3.7) 
is valid. Still, a phenomenologically relevant range of parameters could be studied; the 



obtained scaling is different from eq. (3.7), but it is still an open question whether the latter 
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is valid or receives corrections due to the failure of the one-dimensional triangular barrier 
approximation. 

Regarding the decays from the ISS vacuum to the vacua involving spectator fields, the 
situation seems less clear. From the scale ratios, only will be relevant, since it influences 
the tunneling along the fields X and Y , whose VEVs are controlled by /i^. However, it is 
hard to guess how dominant the /x^-dependent pieces in the potential are in comparison with 
other terms, and if at all, the decay rate is expected to increase for growing /x^, since the 
difference between the VEVs in the ISS and spectator vacua grows with Concerning A 
and ka, although they control the interaction with the spectator fields, typically one has 



K\ ^ 1, which implies that the VEVs of the spectator fields will be very small -see eq. (2.9); 
thus the contribution of their kinetic energy to the bounce action will be sub dominant, and 
their profile in the bounce configuration will be approximately that which minimizes the 
potential for the rest of the fields at any given point in the bounce trajectory. The value 
of the potential at these minima won't be affected by changes of A, ka, and so they will 
have a weak influence in the final value for the bounce action. Thus, one expects that the 
latter will be mainly determined by the remaining dimensionless parameter which as said 
before cannot be made parametrically small for phenomenological reasons. Therefore, a 
more detailed numerical evaluation is needed in order to elucidate whether sufficiently long 
lifetimes can be achieved for phenomenologically acceptable values of the parameters. 

The following sections present the results of the numerical evaluations of the bounce 
configurations and their actions. Given the preceding discussion, the computations are done 
for different values of ka, and h, keeping fixed A^c = 5 as well as A^c = 1 and A = 1. 

3.3 Numerical results: decays towards susy vacua with no R-symmetry break- 
ing 

This section presents the results for the numerical evaluation of the bounce action associated 
with transitions from the ISS to the susy vacua generated by nonperturbative interactions, in 
the case /i,^ = 0. Both vacua differ in the VEVs of the fields X, F , x, in both of them one 
can choose (x) = (x); which will be respected by the evolution under the equations of motion 
of the bounce. Turning on other fields apart from the above will increase the potential energy 
by turning on additional F-terms; thus the bounce configuration will only involve nontrivial 
profiles for three fields: X, Y and x = X- The parameters in the superpotential are taken 
as real, which allows to only consider tunneling along the real parts of the above fields; 
furthermore, X is treated as a diagonal field, X = xl^. In order to solve for the bounce 



configuration with the technique of ^3.1, the fields have to be rescaled so that their kinetic 



terms are canonical as in eq. (3.2); thus, the following redefinitions are introduced. 



x=^x', Y = -Y', x = X=-x'- (3.8) 



The bounce action was computed for different values of the parameter ka, which according 
to the discussion in the previous section should largely determine the decay rate. The result 
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is shown in Fig. [TJ which confirms the parametric dependence in ka, with exponent equal 
to 2.6. This confirms the analytical result of eq. (3.6), which for the choice of parameters 



here yields an exponent of 2.4. The results show that the normalization is such that even for 
moderate values of ka the lifetimes are much bigger than the age of the Universe, Sb 3> 400. 
One of the resulting multifield bounce configurations is shown on Fig. [2j 




400000 



300000 - 



200000 - 



Figure 1: Bounce action corresponding to the decay towards the susy vacua generated by 
nonperturbative interactions, for different values of kx, = 5, Nc = 1, = 1,k^ = 0. The 
solid line represents the function 5*;, 



3324-6 + 2068. 



3.4 Numerical results: decays towards susy vacua with R-symmetry breaking 

By similar arguments to those employed in the previous section, one can conclude that the 
bounce configuration will again include only nontrivial profiles for the fields x = X^^y ^• 



Again, the rescalings of eq. (3.8) are used, and the bounce action is computed for different 
values of the parameter k^, keeping = 5, iVc = l,h = 1, as well as choosing a very 
high value of ka = 3.85 ■ 10^ in order to make the nonperturbative interactions irrelevant. 
The results are shown in Fig. [3} which shows that the bounce action increases very rapidly 
for small k^, following a power law with exponent approximately equal to 8, which differs 



from the estimate of eq. (3.7). This estimate is only valid for k,^ <^ 1, which could not be 
probed due to numerical fluctuations; hence the discrepancy should not be too surprising. 
An example of a bounce configuration is shown on Fig. Ill Though the scaling with differs 



from that of eq. (3.7), qualitatively the bounce configuration behaves as described in ref. [2]: 
x' and Y' only have small variations and the tunneling proceeds mainly in the x' direction. It 
should be noted that the enhanced growth of the bounce action with diminishing means 
that the lifetime of the ISS vacuum will be long enough for typical phenomenologically 
acceptable models. In these, /z^ fixes the gaugino masses while fi determines the masses of 
the messenger fields, which gives /x<^ ~ ITeV, /i ~ 200TeV and hence ~ 0.005, which is 
the region studied here, which yields Sb ^ 400. Thus one expects no phenomenologically 
relevant constraints coming from demanding sufficiently suppressed decays. 



13 



x' 




Figure 2: Bounce configuration and potential corresponding to the decay towards the susy 
vacua generated by nonperturbative interactions, for Nc = 5, Nc = ^, h = 1, = 0, kx = 15. 



3.5 Numerical results: decays towards vacua with nonzero VEVs for spectator 
fields 

As argued in §3.2[ in the presence of spectator fields the decays of the ISS vacuum are 
expected to proceed dominantly towards the vacua involving nonzero VEVs for diagonal 
spectators. This section presents results concerning decays towards a vacuum in which all 
the spectator fields coupling to the diagonal of the meson field X get VEVs, for Nc = I and 
Nc = 5. The fields that get different VEVs in the two vacua are the diagonal components 
of X and the spectator S, as well as Y and p^.p^. Again, boundary conditions motivate us 
to consider x = x as well as p = p along the bounce trajectory. Once more we take X as 
diagonal, X = xl^, and the spectator field matrix is taken as S" = diag(s, s, s, s, — 4s). To 
get canonical kinetic terms, the fields have to be rescaled as 

X = x', Y = - Y', y = y = - y', = = -p', s = —= s' . 
yio 2 ' ^ ^ 2 2 ^40 



Thus, the bounce configurations involve 5 fields, x' ,Y' ^x' ,r' and s' . The bounce action was 
computed for varying h as well as for varying k^, keeping the rest of parameters fixed; the 
results are shown in Figs. [5] and |6} respectively. Fig. [Tjshows one of the bounce configurations. 
From Figs. [5] and [6] it is clear that one can obtain long-lived enough vacua with 5*;, > 400 for 
order one values of h, which confirms the result quoted in ref. |6], which were obtained from 
a simpler calculation which ignored the tunneling along the X and Y directions. The more 
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Figure 3: Bounce action corresponding to the decay towards the susy vacua generated by 
R-symmetry breaking interactions, for different values of k^, Nc = 5,Nc = l,h = 1, ka = 
3.85 ■ 10^ The sohd hue represents the function St = 7.25 ■ IQ-^^/t/-^^ - 2.98 ■ 10^. 



refined calculation presented here shows that the effect of turns out to be very relevant, 
as shown in Fig. [6] The bounce action scales with in a predominantly quadratic way. The 
results suggest then the dependence given in eq. (1.1), valid for the choice of parameters of 
this section. 
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A One-loop Coleman- Weinberg potential 

This paper presents computations of tunneling configurations using as potential the tree- 
level one plus the one-loop contribution to the effective potential, which is obtained as the 
Coleman- Weinberg potential Vcw evaluated for arbitrary background values of the fields. 
Using the supersymmetry preserving DR scheme, the only contributions to Vcw come from 
integrating out the fields whose spectrum is not supersymmetric at tree-level. Since we 
are interested in the effective potential along the bounce trajectories interpolating between 
different vacua, we should take into account the fields that get a nonsupersymmetric spectrum 
not only in one of the vacua but also along the said trajectory. 

In the ISS vacuum, the fields with a nonsupersymmetric spectrum are p, p, Z and Z, 
while in the vacua associated with the vector spectators, these are a combination of x, 
some of the components of p, p and off-diagonal components of X. Gauge fields should also 



15 



2.00 h 



1.99 
1.98 
1.97 
1.96 ^ 



20 40 



80 100 



20 40/ 60 80 100 



Figure 4: Bounce configuration and potential corresponding to the decay towards the susy 
vacua generated by R-symmetry breaking interactions, for A^^ = 5, iY^ = 1, /?. = = 
0.005, = 3.85 ■ 10^. 



be taken into account: a crucial ingredient in the stabilization of all field directions in the ISS 
vacua is the gauging of part of the global symmetry, which generates D-term contributions 
that can break the degeneracy in the vector supermultiplet. This does not happen in the 
metastable vacua described in ^but it may in the trajectories interpolating between them, 
so that new contributions to the one-loop potential might arise, in particular those yielding 
the mass terms that stabilize the pseudo-Goldstone modes of the ISS vacuum. These were 
modeled in ref. |6j as a term of the form 

VD-^\p + p*\'\x + X*\^ 

where g is the coupling of the gauged symmetry -note how these contributions vanish in 
both the ISS and spectator vacua, with p = p = and X = X = 0, respectively. In the 
computations presented here these terms are not put by hand but will be encoded in the 
gauge field contributions to the Coleman- Weinberg potential; for simplicity a U(l) gauge 
symmetry is considered, under which only the fields p, p are charged with charge -|-1 and — 1, 
respectively. 

In order to make the computation feasible by having analytic expressions for the mass 
matrices, it was decided to ignore the contributions from massive components of X to the 
Coleman- Weinberg potential. This is only relevant when considering decays from the ISS 
vacuum to spectator vacua, but it is a good approximation due to the fact that the bounce 
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Figure 5: Bounce action corresponding to the decay towards the vacua with nonzero spectator 
VEVs, for different values of /i, = 5,Nc = 1, ^^a = 38.5, k,^ = 7.7 ■ 10""^. The sohd hne 
represents the function 160.728h-^-^^^^^ - 9.79449. 



configurations exit the barrier close to the ISS vacuum, in which the field components of the 
superfield X have a supersymmetric tree-level mass matrix. 

The one-loop Coleman- Weinberg potential at a scale Q in the DR scheme can be written 

as 

n 

where Sn = 1, 1/2, for vector, fermions and scalars, respectively, and m„ are the correspond- 
ing field-dependent mass eigenvalues in a proper normalization of the fields. As justified in 
the discussion above, we include the following degrees of freedom: 

• Vector: U(l) gauge field, 

• Fermions: U(l) gaugino, fermionic components of the chiral fields p, p, x, Xi ^i ^i 

• Scalars: Scalar components of the chiral fields p, p, Xi Xi ^? ^■ 



The mass eigenvalues in eq. (A.l) should be obtained for arbitrary background values of 
the fields with nontrivial tunneling profiles; by looking at the tree-level potential, it can be 
justified that these will be the fields taking nonzero VEVs in either of the vacua involved 
in the decay, since considering nonzero VEVs for additional fields will always increase the 
potential energy and drive the field configuration away from the sought-for critical point 
of the action. Thus, the mass matrices of the fields mentioned above will depend on the 
background values of the fields Xi Xi S, X, Y and p, p. Given the boundary conditions for 
the semiclassical considerations, one can restrict to configurations with p = p, X = X- Fur- 
thermore, we will consider diagonal backgrounds for X and S, and restrict to Nc = 5, so 
that X = xls, S = diag(s, s, s, s, —4s). In order to write down the mass matrices, since we 
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Figure 6: Bounce action corresponding to the decay towards the vacua with nonzero spectator 
VEVs, for different values of /t,^, Nc = 5, Nc = ^, h = 0.6, ka = 38.5. The sohd hne represents 
the function 8.9719 x 10^4 + 153663.^0 + 321.42. 



consider vacua in which one out of the components of p, p gets a VEV, we separate the 
p, p fields and the Z,Z with which they mix in two groups, p = (pi,p2), where pi is the 
component getting a VEV in the spectator vacua, and similarly for p, Z, Z. We then define 
the following multiplets: 

Vector : b^, 

Fermion : ^ = (x, pi, Zi, x, Pi, ^i. A), ^' = (pa, Z2, P2, ^2), 
Scalars : S*! = Re(x, Pi, ^1, X, Pi, ^1), 5'2 = Im(x, Pi, ^1, X, Pi, ^1), 
S[ = Re(p2, Z2, P2, Z2), S2 = Im(p2, Z2, P2, Z2), 

where 6^ and A are the U(l) gauge field and gaugino. The mass matrices in the normalization 
consistent with eq. ( A.l ) are, for arbitrary background values of the fields X,Y, pi = pi,x = 
X, s, X, and for Nc = 1, 

j=l,2 
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hx hx 

hx 

hx hx 

hx h'^fJ'4, 
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{Ml,),, = {Ml,),, = 

mil = "^44 = h^p\ + hY"^, mi2 = 17145 = 

mi3 = m46 = h^l^^Pi + /i^Pii^, mi4 = h'^{x^ - l^^) + h'^P^X, 

mi5 = m24 = 0, mi6 = ^34 = h'^xpi, 

= /i^(Pi - A«^) + /i^/^^a; - Sg-^pi - 4AA/is, m26 = TO35 = /i^l^x, 

(^1,2)^ = (^l,2)ii = riij, 

nil = ?T'44 = h'^Pi + hY"^, ni2 = 7145 = 0, 

7^13 = 7146 = /i^/^^Pi + h'^PiY, nu = ^^(x^ - A*^) + ^^A*.)!.^, 

7115 = "-24 = 0, 7li6 = 7l34 = h'^Xpi, 

7l22 = 7l55 = h'^X^ + h'^X^, 1^,23 = 7l56 = /i^A<^X + h'^X^, 

7125 = ^^(P? - 1^'^) + h^l^,pX - AXAhs, 7l26 = ^^35 = h^^X, 

nss = nee = h^pl + hi^pj + /^^x^ "-36 = "-66 = 0, 
M'li = M'l^ = 

/i^X^ + h^x^ h^l^ibX + /i^xa^ —h?p? + /i^/x</,a; + /lAAs h^Yx 

h^P4,X + ^'^X^^ ^"^A^^ + h'^X^ h'^^X 

-/i^/x^ + h^ii^x + /iAAs /i^i^X h'^X^ + ^^'"^^ h^l^4>X + /i^xa^ 

/i^Fx h^lJ'<t>X + /i^xa; ^^A''! + ^^X^ 
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